Social environment influences microbiota and potentially pathogenic bacterial communities on the skin of developing birds

Background Animal bacterial symbionts are established early in life, either through vertical transmission and/or by horizontal transmission from both the physical and the social environment, such as direct contact with con- or heterospecifics. The social environment particularly can influence the acquisition of both mutualistic and pathogenic bacteria, with consequences for the stability of symbiotic communities. However, segregating the effects of the shared physical environment from those of the social interactions is challenging, limiting our current knowledge on the role of the social environment in structuring bacterial communities in wild animals. Here, we take advantage of the avian brood-parasite system of Eurasian magpies (Pica pica) and great spotted cuckoos (Clamator glandarius) to explore how the interspecific social environment (magpie nestlings developing with or without heterospecifics) affects bacterial communities on uropygial gland skin. Results We demonstrated interspecific differences in bacterial community compositions in members of the two species when growing up in monospecific nests. However, the bacterial community of magpies in heterospecific nests was richer, more diverse, and more similar to their cuckoo nest-mates than when growing up in monospecific nests. These patterns were alike for the subset of microbes that could be considered core, but when looking at the subset of potentially pathogenic bacterial genera, cuckoo presence reduced the relative abundance of potentially pathogenic bacterial genera on magpies. Conclusions Our findings highlight the role of social interactions in shaping the assembly of the avian skin bacterial communities during the nestling period, as exemplified in a brood parasite—host system. Supplementary Information The online version contains supplementary material available at 10.1186/s42523-024-00327-2.

The social environment (i.e., environmental characteristics of interacting individuals) has been suggested to promote similarities in bacterial communities between interacting animals [28][29][30][31].Thus, it should be important in driving similarities in microbiome-derived physiological and behavioural traits of hosts, as well as explaining susceptibility to parasitism [1,35].Most evidence for effects of the social environment on the microbiome comes from experimental approaches in a few captive animal models, or from correlational studies in gut and skin microbiota of humans [31,[36][37][38], non-human primates [28,29,32,[39][40][41], other mammals [42,43], birds [44][45][46][47], amphibians [48], and arthropods [49,50].Despite evidence supporting associations between social interactions and the host microbiota, we are just starting to understand how social interactions structure host microbial communities.This is partly due to the confounding effects of sharing environments in the absence of social interactions and the role of host genetics [32].Within social groups, individuals are likely to share early life environmental conditions, physiological stress, similar resources, diet, and/or genetic relatedness (family groups).Similarities in microbial communities among individuals would be predicted in such cases even in the absence of social interactions.Thus, to understand the role of social interactions, it is essential to disentangle these effects, and studying unrelated individuals interacting under identical environmental conditions may help achieve this [51][52][53][54].
In the present study, we take advantage of the brood parasite -host system formed by great spotted cuckoos (Clamator glandarius) (hereafter cuckoos) and Eurasian magpies (Pica pica) (hereafter magpies).Cuckoos are obligated brood parasites that in Europe mainly lay their eggs in magpie nests, where magpie adults incubate the eggs and take care of the cuckoo chicks during the nestling and fledgling periods [55].Adult cuckoos thus do not have contacts with their own nestlings, restricting microbial transfers from parents to offspring to the prelaying phase.When the eggs hatch and parasitic nestlings do not outcompete host nestlings, host adults rear their own nestlings along with parasite nestlings [55].In these cases, the skin and feathers of parasitic and host nestlings are in close contact, which allows exploring similarities in bacterial communities of natural and foster siblings that cannot be explained by relatedness.During the nesting phase, nestlings of both species share similar environmental conditions, including those related to parental care.However, magpie nestlings that develop together with cuckoos differ in the social environment from those that grow up in monospecific nests.Thus, consistent differences in the microbiota of host nestlings that do or do not share nests with cuckoos can be interpreted as the results of social interactions with heterospecifics.
Capitalizing on this natural system, we conducted a cross-fostering experiment, where we manipulated the heterospecific social environments to disentangle the effect of physical and social environments on assembly processes of uropygial gland skin microbiomes (Fig. 1).The cross-fostering approach allowed us to avoid possible biases due to cuckoos choosing nests of particular environmental characteristics [56,57], and to maximize the number of magpie nests with individuals of both species.Our experimental design allows testing effects of exposing individuals of the same social group (here members of the same species: magpies) to members of a different social group (cuckoos) in the same environment (magpie nests and parents).We characterized the bacterial community of the uropygial gland skin of magpies and cuckoos developing in con-or heterospecific broods via amplicon sequencing of the bacterial 16S rRNA gene.We focused on the uropygial gland skin, as this gland produces a secretion that birds spread onto their feathers and skin while preening [58], and from which several bacterial strains have been isolated [59][60][61][62][63][64][65][66][67].The uropygial gland secretion has a species-specific chemical composition [58] which, in contact with skin and other body parts, may act as a substrate that improves the establishment and growth of species-specific bacterial communities [68,69].Therefore, we first hypothesised that the skin microbiota of chicks from magpie nests with only magpies or cuckoos (monospecific nests) would vary interspecifically, because of the impact of host intrinsic characteristics on the microbiota (Fig. 1A).Secondly, assuming transmission of microbes between nest-mates in a shared social and physical environment, we expected that chicks growing in heterospecific nests would show reduced interspecific differences (Fig. 1B).Third, given the expected effects of social transmission, we hypothesised that microbial communities would differ between magpie chicks that did or did not grow up with cuckoos (Fig. 1C).We explored these hypotheses in the skin bacterial community as a whole, but also in subsets constituting the core microbiota and a set of potential pathogens.

Study area and fieldwork
Fieldwork was performed during the breeding season of 2018 in the Hoya de Guadix (37°15′N; 3°01′W), a semiarid high-altitude plateau in southern Spain, where a magpie population frequently parasitized by cuckoos is established [70,71].The vegetation is sparse, with disperse holm-oak trees (Quercus rotundifolia), groves of almond trees (Prunus dulcis) and pines (Pinus halepensis), where magpies usually build their nests.
Since mid-March, we intensively searched for magpie nests, which allowed us to infer the start of laying (hereafter laying date) and, thus, the expected hatching dates of cuckoo and magpie eggs.Cuckoos might choose to parasitize host nests with characteristics that maximize their reproductive success [56], which might result in biased samples when exclusively using natural parasitized nests, and, thus, highlight the importance of adopting an experimental approach.Moreover, magpie nestlings usually starve in naturally parasitic nests because cuckoo eggs hatch four or more days ahead of those of magpie, which are outcompeted by older cuckoo chicks [72].To avoid possible bias of naturally parasitized nests, and to maximize the probability that magpies and cuckoo nestlings develop together in the same magpie nests, we cross-fostered cuckoo eggs between magpie nests whenever possible, synchronizing their expected hatching date of cuckoo eggs to the same or one or two days later than that of the magpie eggs.This approach allowed us to maximize the number of heterospecific nests where cuckoo and magpie nestlings developed together [see 19 for a similar experimental approach].Magpie monospecific nests were simply non-parasitized nests, whereas cuckoo monospecific nests were parasitized magpie nests where natural death of magpie nestlings occurred (Fig. 1).Cuckoo and magpie nestlings develop at different rates [73] and, thus, to sample them at similar developmental stages it was necessary to visit the nests twice.However, we only manipulate the species needed for processing at each visit (either magpie or cuckoo nestlings), leaving the nestlings of the other species in the nest.We collected microbiome samples from 56 nestlings when cuckoo and magpies were, approximately, 15 and 17 days old respectively.Briefly, we sampled skin microbiota of nestlings by rubbing the surface skin of the uropygial gland, including the opening, with a sterile cotton swab (APTACA, ref. 2160, Canelli, Italy) wetted in sterile Phosphate Buffer Saline (PBS, 0.2 M).The swab with the bacterial sample was kept in a sterile microfuge vial with 1 mL of sterile PBS and stored at -18 °C until DNA extraction.At the time of sampling, we also measured tarsus length with a digital calliper (accuracy 0.01 mm), wing length with a ruler (accuracy 1 mm), body mass with a digital scale (accuracy 0.01 g), and gland dimensions (length, width and height) with a digital calliper (accuracy 0.01 mm) as described previously [74].

DNA extraction and amplicon sequencing
DNA from the sampled bacterial communities was extracted using the FavorPrep ™ Blood Genomic DNA Extraction Mini Kit (Favorgen Biotech Crop., Taipei, Taiwan), following this protocol: first, the sample was sonicated for 2 min at 120 Hz to release the bacterial Fig. 1 Experimental approach scheme, including the types of magpie nests in the study as well as sample sizes of magpies and cuckoos in each type of nests.The predicted similarities ( =) and differences ( ≠) in microbial composition between nestlings are indicated with green labels cells from the swab.The swab was then removed, and the PBS with the bacteria was centrifuged at 13,000 rpm for 5 min.The supernatant was discarded, and 180 µl of TES (25 mM Tris-HCl, pH 8, 10 mM EDTA and 10% sucrose), 10 mg/ml of lysozyme and 10 mg/ml of RNase were added to the pellet.Subsequent steps were performed according to the FavorPrep ™ protocol.From this extraction, 5 µl were used to perform a PCR reaction to verify the presence of bacterial DNA in the uropygial skin.The PCRs were conducted using the primers B969F (ACGCGHNRAA CCT TACC) and BA1406R (ACG GGC RGT GWG TRCAA) [75].The PCR products were visualized on a 1% agarose gel with electrophoresis.Libraries for paired-end Illumina sequencing were constructed in two steps following Caporaso approach [76] with the bacterial V6-V8 region of the 16S rRNA using the same pair of primers specified above.These primers maximize the amplification of bacteria and reduce non-specific eukaryotic amplifications [77].Library construction and paired-end sequencing (2 × 300) in the MiSeq (Illumina) platform was carried out at the Institute of Parasitology and Biomedicine "López-Neyra" facilities (IPBLN, Granada, Spain).

Bioinformatic amplicon data processing
We first processed the amplicon sequences in QIIME2 v2020.6 [78], using default parameters unless stated otherwise.Primer trimming and sequence quality filtering were performed using DADA2 [79], and all sequences were clustered into ASVs (Amplicon Sequence Variants) at 100% similarity and assigned to taxonomy using the Silva 138 database [80].Due to the primers' specificity for bacteria, non-bacterial sequences, and sequences identified as mitochondrial or chloroplast, were removed from the ASV table.Contaminant sequences were identified from field (open swabs without sample) and laboratory (extraction and sequencing blanks) negative controls with the "Decontam" package in R [81,82] using the prevalence method and a threshold of 0.4.Sequences were aligned and a rooted bacterial phylogeny was generated using the method align-to-tree-mafft-fasttree in QIIME2.One sample was filtered out due to markedly low reads (2,500).The ASV table was rarefied to the minimum sampling depth (14,877 sequences) using the method rarefy_ even_depth in the "phyloseq" package [83].
We also identified potential avian pathogens and the core microbiome for each species among the detected ASVs.For the potential avian pathogens, we first ran the FAPROTAX script in python [84], which converts prokaryotic abundance tables (ASV tables) into putative functional abundance profiles.The ASVs considered by FAPROTAX as animal pathogens were searched in the literature to certify avian pathogenicity.Besides, we also used the Pathogen Host Interaction database (PHI-base) [85] and the review published by Benskin et al. in 2009 [86] to search for genera that includes potential known pathogenic bacteria of birds.We use those datasets to build a new ASV table that included potential pathogenic ASVs belonging to genus with available information (Additional file 1).We also calculated the core microbiome using a relative abundance of 0.0001% in at least 50% of the samples in the "phylosmith" package [87] in R v4.0.2 [82].We did so separately for each species and type of social environment considered (i.e., only magpies, only cuckoos, magpies that grew with cuckoos and cuckoos that developed together with magpies).Then, we created a subset of the ASV table pruning out taxa that did not belong to the core microbiome for each species.
Factors expected to influence alpha and beta diversity indexes were respectively explored in mixed model ANOVAs and PERMANOVAs.The effects of species identity (hereafter, ID) were explored with information from nests where only cuckoo or only magpie nestlings developed (monospecific nests).The models included species ID as fixed factor and the nest ID (nested within, species ID) as the random factor.The effects of species ID were also explored in nests where magpie and cuckoo nestlings develop together (heterospecific nests), but in this case, the statistical model included species ID as the fixed factor and nest ID and the interaction of nest ID with species ID as random factors.The effect of social environment was analysed by comparing magpie nestlings that grew in monospecific nests with those that developed together with cuckoos in heterospecific nests.These models included the social environment (mono-or heterospecific magpie nests) as the fixed factor and nest ID (nested within social environment) as the random factor.Brood size did not significantly explain alpha diversity indexes (Additional file 2) and, thus, was not included as covariable in the statistical models.
We tested which bacterial genera had significant differential abundances among the four types of social environments.We did this by using the trans_diff function from "microeco" package [90] in R v4.0.2 [82] with the Wilcoxon Rank Sum method and False Discovery Rate (FDR) adjusted p-values.We conducted differential abundance analyses with the whole ASV table.Finally, we used betadisper function in the "vegan" package [91] using spatial median and adjusted biases to analyse the homogeneity of variances among magpies sharing and not sharing nests with cuckoos.The effects of species ID and social environment on alpha and beta diversities of subsets that included potentially pathogenic bacteria, or the core microbiome, were explored in statistical models identical to those described above.ANOVAs were conducted in STATIS-TICA v.12 [92], while PERMANOVAs were performed with Primer7 v.7.0.17 (PRIMER-e).

Microbial diversity
Considering the whole set of ASVs, alpha diversity indexes of magpie and cuckoo samples did not differ significantly (Table 2).That was the case independently of comparing samples from mono-or heterospecific nests (Table 2).However, alpha diversity of magpie samples from heterospecific nests was significantly higher than that of magpie samples from monospecific nests (Table 2, Fig. 2B).Interestingly, magpies growing up in heterospecific nests shared more ASVs with cuckoos than with magpies or cuckoos from monospecific nests (Fig. 2C).
When considering the beta diversity of cuckoo and magpie samples from monospecific nests, the bacterial communities differed significantly between the two species [except for weighted UniFrac distances (Table 3)], and their 95% confidence interval ellipses in PCoA plots hardly overlap (Fig. 3A).This effect disappeared when comparing samples of magpies and cuckoo nestlings that were raised in the same nest (Table 3), as revealed by overlapping points and 95% confidence intervals in Fig. 3A.Furthermore, regardless of the distance matrix used, the social environment influences the composition of the bacterial community of magpies (Table 3, Fig. 3).However, the individual variation in magpie microbiomes was not associated with social environment (Fig. 3C) when considering Bray-Curtis, Jaccard, or weighted UniFrac distance matrixes (betadisper test; F 1,38 < 2.64, p > 0.109), but it was associated with social environment when considering Unweighted UniFrac (betadisper test; F 1,38 = 23.21,p < 0.001).
Thirty bacterial genera were significantly differentially abundant between groups (Additional file 5).Thirteen, 17 and 21 genera were respectively more abundant in cuckoos from monospecific nests, cuckoos from heterospecific nests, and magpies from heterospecific nests than in magpies from monospecific nests.Magpies from monospecific nests had higher relative abundances of Enterococcus than magpies from heterospecific nests, and higher relative abundances of Escherichia -Shigella and Pseudomonas than magpies and cuckoos from heterospecific nests (Fig. 3D, Additional file 5).Interestingly, magpies in heterospecific nests did not present any differentially abundant genera compared with cuckoos, either from mono-or heterospecific nests.

Core microbiome
Core microbes with relative abundance higher than 0.0001% and that appeared in 50% of the samples for a species comprised 232 ASVs from 75 genera in six phyla (Additional file 6).For these cores, we detected interspecific differences in alpha and beta diversity between nestlings of the two species from monospecific nests (Tables 2 and 3).Furthermore, and in agreement with the results for all ASVs, bacterial diversity and community composition of birds in heterospecific nests did not differ between species (Tables 2 and 3).In addition, differences in social environment for magpie nestlings were also apparent when exploring the core microbiome (Tables 2  and 3).

Diversity and composition of potentially pathogenic ASVs
Alpha diversity of potential pathogenic taxa did not differ significantly among magpies and cuckoos when considering either mono-or heterospecific nests (Table 2).Similarly, sharing the nest with heterospecifics did not affect the alpha diversity indexes of magpie samples (Table 2).When looking at beta diversity of communities of potential pathogens, we found statistically significant interspecific differences when comparing monospecific nests and considering Bray-Curtis or Jaccard distance matrices (Table 3).However, this effect disappeared when comparing cuckoo and magpie nestlings that grew up within the same nests (Table 3).
Moreover, beta diversity of potential pathogenic bacteria of magpies differed according to social environment when considering Bray-Curtis and Jaccard distance matrixes (Table 3).Furthermore, Pseudomonas, Escherichia-Shigella and Enterococcus were significantly more abundant in magpies that grew up in monospecific nests (Fig. 3D).

Discussion
By capitalising on the natural associations between magpie hosts and cuckoo brood parasites, we document the influence of social environment on microbial assemblages of the uropygial gland skin of developing magpie and parasitic cuckoo chicks.As expected, we detected species-specific diversity and composition in both complete microbiomes and microbial cores, underlining the role of species-specific factors in shaping microbial assemblages.However, and as predicted, these interspecific differences disappeared when species cohabitated the same nests, implying either social transmission of symbionts among nestlings or transmission of microbes via feeding adults and/or the shared physical environment [c.f., 93].In heterospecific nests, we observed a change in the magpie microbiota to resemble that of their heterospecific nestmates that most likely imply horizontal transmission of microbes between cuckoos and magpies via their altered social environment.Lastly, we observed a negative impact of social environment on certain candidate pathogenic bacteria in magpies, suggesting potential positive effects of the presence of cuckoos on magpie skin microbiomes.
We found interspecific differences in the uropygial gland skin microbiomes in monospecific nests when using any of the distance matrixes, with the exception of Weighted UniFrac.UniFrac distances are controlled for phylogenetic association of the considered ASVs, and the lack of differences may suggest that the detected interspecific differences for other distance matrices are driven by the relative abundance of phylogenetically closely related ASVs.These interspecific differences cannot only be Table 2 Results from mixed model ANOVAs exploring the effects of species ID in either, mono-or heterospecific nests, as well as the effect on social environment on the alpha diversity indexes (Shannon's diversity index and Faith's phylogenetic diversity (PD)) on bacterial communities of the uropygial gland skin of magpies and great spotted cuckoos.We analysed the effects of species ID and social environment on diversity of the whole bacterial community, the core microbiome and the subset of potentially pathogenic ASVs.Results with associated p-value lower than 0.05 are shown in bold.Fixed (f ) and random (rnd) factors are indicated   accounted by vertical transmission of microbes from the biological or foster parents.This is particularly true for cuckoos, as vertical transmission of microbes is restricted to the pre-egg laying period [94,95].These differences may thus result from intrinsic factors, such as for example species-specific chemical composition of the uropygial gland secretion [58,68,69].Even during the nestling stage, birds preen their feathers and skin with uropygial secretion [96], and thus the antimicrobial properties of these secretions [97,98] may prevent specific bacteria from establishing [99][100][101], and stimulate the growth of other microbial taxa [69].Consequently, it is likely that the particularities of the uropygial secretion of magpies and cuckoos promote a species-specific selective environment that favours certain microbes to grow on the skin, a possibility worth exploring in the future by testing for promoting or inhibitory effects of the uropygial secretion on the bacterial strains.Another non-exclusive possibility explaining interspecific differences is the potential transfer of species-specific fecal microbes [19] to the skin of the uropygial gland.Cuckoos have defensive faeces with a strong smell to deter predators [102], which might host particular bacterial taxa that could be responsible for the detected interspecific differences.Some of these chemical-producing bacteria are likely anaerobic bacteria, which would explain the higher prevalence of the anaerobic Bacteroides, Clostridium, Parabacteroides, and Lachnoclostridium in cuckoo samples.
In gut microbial communities of cuckoos and magpies, species specificities in cloacal microbiomes are retained in heterospecific nests [19,21].In contrast, we found that skin microbiomes converged in the host-parasite species pair.This indicates that the magnitude of the effect of social and shared physical environment varies depending on whether the host-associated microbiomes are external or internal [34,103].Despite similarities in the microbial composition (e.g., saliva microbiomes) and diet between cuckoo and magpie nestlings in the same nests [19], species-specific cloacal microbiomes imply internallymaintained digestive tract microbiomes that are resistant to perturbations from the shared environment [21], while skin microbiomes are more susceptible to horizontal transfer.
The evidence for social transmission of microbes stems mainly from the potential transfer of cuckoo-specific microbes to magpies.This unidirectional transfer of microbes from cuckoos to magpies might indicate that the resilience of skin microbiomes to the social environment also varies interspecifically.Given the broodparasite lifestyle, cuckoos may depend more on in ovo vertical transmission of microbes than magpie nestlings [c.f., 104], while having more resistant skin microbiomes to ensure transgenerational transfers of symbionts [c.f., 93].This skewed opportunity for vertical transmission of microbes in cuckoos may alter host-symbiont associations in this species with potential implications for losses or replacements of microbial symbionts across generations.However, the effects of sharing a nest with heterospecifics might only be temporal, as shown by a cross-fostering experiment on captive zebra finch (Taeniopygia guttata) and Bengalese finch (Lonchura striata domestica) nestlings [45].In these cases, effects of foster families on cloacal microbiomes early in the nestling period disappeared in later stages [45].In the case of the uropygial gland, its secretion might not be fully developed during the nestling stage [58] and, thus, the associated microbiome is likely shaped during the second part of the nesting phase [105], explaining why we detected the expected effect of social environment at the late developmental stage of magpie nestlings.However, to fully grasp the breadth of how early life social environment influences long-term associations and generational transfers of skin symbionts, we need to explore the fate of skin microbes over time during an individual's life.
The social transfer of potentially pathogenic bacteria contrasted the patterns for non-pathogenic bacteria, where the relative abundances of potentially pathogenic genera Pseudomonas, Escherichia-Shigella and Enterococcus were significantly higher in magpies in monospecific than heterospecific nests.The lower prevalence of these genera in heterospecific nests could be mediated by the parallel increase in the diversity of the skin microbiota of nestlings.This might be because increased microbial diversity provides increased resistance to pathogen colonization [106,107] and also stimulates the host immune system [11].Alternatively, the properties of the uropygial secretion, or symbiotic defensive bacteria within the gland of cuckoos, may counter potential pathogens.Even if this might have a positive effect of reduced pathogens on host magpies, it would be unlikely to counter the negative fitness effects of brood parasitism [108].Nevertheless, because we did not test for pathogenicity of strains, these inferences are only tentative and, to fully understand the implication of social environment on pathogenic taxa, future research is needed to explore the specific bacterial strains with detrimental effects on birds, and how they are distributed among con-and heterospecific nests.Furthermore, to disentangle the mechanisms underlining the observed negative effect of social environment on potentially pathogenic bacteria, we need to isolate these bacterial taxa and conduct coculture assays with cuckoo and magpie uropygial gland secretions.

Conclusions
Using a natural host-brood parasite system, we were able to separate the effects of genetic relatedness and shared environment on the skin microbiota from those of interacting individuals, elucidating a role of social environment determining the skin microbiomes of wild birds.Our study implies that skin microbiomes are amenable to horizontal transfer of microbes from the social and the nest environment, but that the magnitude and identities of bacterial genera transferred depend on host ecology.Early-life exposure to heterospecific microbes can thus alter wild bird skin microbiomes, which should influence both short and long-term stability of beneficial and antagonistic symbiotic interactions.

Fig. 2 A
Fig. 2 A Microbial composition at the phylum and genus levels of uropygial gland skin of great spotted cuckoos and Eurasian magpies from monospecific (M) or heterospecific (H) nests.B Least square means (± 95% CI) of alpha diversity indexes (Shannon's diversity index and Faith's phylogenetic distance (PD)) estimated for magpies and cuckoos from monospecific (M) or heterospecific (H) nests.C Venn diagram showing the number of shared ASVs between different treatment groups

Fig. 3
Fig. 3 PCoA plots show the microbiome composition based on Bray-Curtis distance matrixes of magpies and great spotted cuckoos from monospecific (M) or heterospecific (H) nests using (A) the whole set of ASVs and (B) only potentially pathogenic ASVs.Ellipses are 95% confidence interval.(C) Box plots showing group dispersions with Bray-Curtis distance matrix (betadisper analyses) of magpies from monospecific (M) or heterospecific (H) nests.(D) Differential abundance of genera from the whole dataset of ASVs differing significantly among species (magpies or great spotted cuckoos) and/or social environments (nestlings from monospecific (M) or heterospecific (H) nests).The list indicates bacterial genera derived from Random Forest analyses with alpha = 0.01.Underlined in red are the potentially pathogenic genera

Table 1
Sample sizes (number of nests and nestlings sampled) for each of the four combinations of species (magpies or great spotted (GS) cuckoos and social environment (monospecific or heterospecific)

Table 3
Results from PERMANOVAs analysing beta diversity matrixes of the whole bacterial community, the core microbiome or potentially pathogenic ASVs subsets as dependent variables.The analyses explored the effect of species ID (taking into account mono-and heterospecific nests), and the effect of social environment on the microbial communities of magpie and great spotted cuckoo nestlings.The number of permutations was set to 9999.Results with associated p-value lower than 0.05 are shown in bold.